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ABSTRACT 

The reasons for the anomalously high thermal stability of cubane CgHg and the mech- 
anisms of its decomposition are studied by numerically simulating the dynamics of this 
metastable cluster at T = 1050 - 2000 K using a tight-binding potential. The decompo- 
sition activation energy is found from the temperature dependence of the cubane lifetime 
obtained from the numerical experiment; this energy is fairly high, Ea = 1.9 ± 0.1 eV. 
The decomposition products dbTQ^ clS Sb rule, either CgHe and C2H2 molecules or the isomer 
CgHg with a lower energy. 
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Cubane CgHg (Fig. 1) discovered in 1964 [1] is of great interest from both the funda- 
mental and practical standpoints. In this cluster, the carbon atoms are in cube vertices; 
i.e., the angles between the C-C covalent bonds are QO'^ rather than 109.5° as in carbon 
compounds with a tetrahedral atomic arrangement and sp^ hybridization of the atomic 
orbitals. Such large bending of the C-C-C bonds is energetically unfavorable. However, 
hydrogen atoms arranged on the main diagonals of the cube stabilize this atomic config- 
uration (Fig. 1), corresponding to a local rather than global minimum of the potential 
energy as a function of the atomic coordinates. Although cubane is a metastable cluster, 
its high stability is demonstrated by the experimental fact that the cubane molecules not 
only retain their structure at temperatures significantly higher than room temperature 
but also can form a molecular crystal, namely, solid cubane s-CgHg, with a melting tem- 
perature near 400 K [2]. The formation heat of cubane is relatively large, 6.5 eV/CsHg 
[3] (such energy is released, e.g., when solid cubane transforms into graphite layers and 
H2 molecules) . High energy content of cubane makes it a promising material for fuel ele- 
ments, and the possibility of replacing the hydrogen atoms with various functional groups 
(such as CH3 in methylcubane [4]) opens the way to the synthesis of new compounds with 
unique properties. 

A wide application of cubane is hampered by the absence of cheap methods for its mass 
production [5]. In our opinion, detailed studies of the mechanisms and products of cubane 
decomposition may suggest a direction of searching for new methods of its production. 

In this case, it is interesting to consider a possible reversal of the chemical reaction 
(e.g., under heating in the presence of appropriate catalysts lowering the barrier of the 
reverse reaction). A certain analogy may be drawn with fuUerene Ceo: its decomposition 
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is preceded by a series of Stone- Wales transformations, which result in the formation of 
"surface" defects and, eventually, a C2 dimer is separated; on the other hand, annealing 
of these defects leads (through the same transformations) to the formation of a fuUerene 
from a strongly distorted spherical cluster Ceo (see [6, 7] and references therein). There 
are only a few experimental studies of cubane decomposition (see, e.g., [8, 9]), and those 
experiments were performed only over narrow temperature [8] and lifetime [9] ranges. 
As for theory, many theoretical studies employed the same schematic potential-energy 
surface of cubane and its isomers [8], while, to the best of our knowledge, the cubane 
dynamics before the instant of transition to another isomer has been studied only at very 
high temperatures and for a very short time (~ 1 ps) corresponding to only several tens 
of cluster oscillation periods [10]. 

The main aim of this work is to numerically simulate the cubane dynamics over a wide 
temperature range and determine the cubane decomposition activation energy, products 
of its decomposition, and the types of isomers forming at the stage of evolution preceding 
the decomposition. We calculated the energies of various atomic configurations within a 
nonorthogonal tight-binding model that was proposed for hydrocarbon compounds in [11] 
and modified in [12] using a criterion of more exact correspondence between the theoretical 
and experimental values of the binding energies and interatomic distances in various 
C„H^ molecules. The model is a reasonable compromise between more rigorous ab initio 
approaches and extremely simplified classical potentials of interatomic interaction. In this 
model, the bond lengths in cubane are calculated to be Ic-c — 1-570 A and Ic-h — 1.082 
A, which are close to experimental values 1.571 and 1.097 A, respectively [13]. The 
calculated binding energy of the atoms in cubane Ef, = [8E{C) +8E{H) — E{C8Hg)]/16 = 
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4.42 eV/atom likewise agrees with the experimental value 4.47 eV/atom [13]. The ratio 
of the energies of the carbon and hydrogen subsystems in the heat-insulatcd cubane 
calculated in this model coincides with the theoretical value [14] . 

We studied the thermal stability of cubane CsHg by the molecular-dynamics method. 
At the initial instant of time, random velocities and displacements are given to each of 
the atoms in such a manner that the momentum and the angular momentum of the whole 
cluster are equal to zero. Then, the forces are calculated acting on the atoms. The classical 
Newton equations of motion are numerically integrated using the velocity Verlet method. 
The time step was to = 2.72 • 10~^^ s. In the course of the simulation, the total energy 
of cubane (the sum of the potential and kinetic energies) remained unchanged, which 
corresponds to a microcanonical ensemble (the system is not in a thermal equilibrium 
with the surroundings [15-18]). In this case, the "dynamic temperature" T is a measure 
of the energy of relative motion of the atoms and is calculated from the formula [19, 20] 
(Ekin) = ^kBT{3n — 6), where {Ekin) is the time-averaged kinetic energy of the cluster, 
ks is the Boltzmann constant, and n = 16 is the number of atoms in cubane. It should be 
noted that the velocity Verlet algorithm is conservative with respect to the momentum 
and the angular momentum [21] and the relative change in the total energy of cubane 
does not exceed 10"*^ for at least 2 • 10^ molecular-dynamic steps, which corresponds to a 
time of ~ 1//S. 

We studied the cubane evolution for 50 various sets of initial velocities and displace- 
ments of the atoms corresponding to temperatures T — 1050 — 2000 K. It turned out that, 
during its decomposition, cubane is most often (in ^ 80% of events) transformed to the 
isomer cyclooctatetraene (COT) (Fig. 2a) with a lower potential energy (a higher binding 
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energy Ej, — 4.82 eV/atom). The cubane decomposition into a benzene CqHq molecule 
{Eh — 4.82 eV/atom) and an acetylene C2H2 molecule {E^ = 4.54 eV/atom) occurs more 
rarely (in 20% of events) (Fig. 2b). We also observed several times the formation of 
styrene and some other isomers of CsHg. Almost without exception, upon decomposition, 
cubane is transformed first into the isomer syn-tricyclooctadiene (STCO) with Ei, — 4.47 
eV/atom (Fig. 3a), which is quickly transformed (in a time of 0.1 - 1 ps) into either COT 
or isomer bicyclooctatriene (BCT) with Ei, — 4.65 eV/atom (Fig. 3b). BCT, in turn, 
transforms into COT or decomposes into benzene and acetylene (BEN + A) molecules. 
We have never observed reverse transitions, such as STCO —>■ cubane, COT — > STCO, 
and BCT — > STCO, whereas sometimes COT is transformed into BCT with subsequent 
decomposition BCT BEN + A. 

As the temperature T decreases from fa 2000 to fa 1000 K, the cubane hfetime r 
increases by six orders of magnitude, from ~ 1 ps to ~ 1/xs (Fig. 4). Since the decompo- 
sition of metastable clusters is an inherently probabilistic process, the lifetime r exhibits 
a dispersion at a given temperature T. Nevertheless, it is seen from Fig. 4 that the results 
of the numerical simulation are described by the common Arrhenius formula 



t-\T) = A-exp 



(1) 



According to this formula, the dependence of ln{T) on 1/T is a straight line, whose slope 
determines the activation energy Ea — (1.9 ± 0.1) eV and its intersection point with the 
ordinate axis determines the frequency factor A = iQ^^-^^^^-^^ g-i. it is remarkable that 
the values of Ea and A agree well with the experimental values Ea — (1.87 ± 0.04) eV 
and A — iq^^-^^^^oaa g-i obtained when studying cubane pyrolysis in a very narrow range 
T = (230 -7- 260) °C [8], which is far apart from the temperature range studied in this 



paper. The small (on the logarithmic scale) difference in the frequency factors is likely 
due to the temperature dependence of A (we note that the value of A for cubane is almost 
four orders of magnitude smaller than that for fuUerene Cgo [7]). 

Extrapolation of the t(T) curve to the range T < 1000 K (which is inaccessible in 
direct numerical calculations because of extremely long simulation time) permits one to 
compare the results of the simulation with the experimental values of r obtained in [9] 
for several temperatures in the range T — 373 973 K. As seen from Fig. 4, here also 
there is agreement between theory and experiment. Thus, Eq. (1) with the values of 
Ea and A found can be used to determine the cubane lifetime (or, in any case, to make 
order-of-magnitude estimates) at both very high and comparatively low temperatures. In 
particular, Eq. (1) gives r ~ 10^^ s at room temperature and r ~ 10® s at the melting 
temperature of solid cubane 400 K (at which, on melting, only weak van der Waals 

bonds between CgHg clusters are broken, whereas the clusters themselves retain their 
structure and the energy stored in them). The lifetime decreases to r ~ 1 s only on 
heating to T Ki 600. Pyrolysis experiments [9] give r 10 ms at T = 573 K and r ^ 2 
ms at T = 673 K (Fig. 4). 

Note that, when analyzing the t(T) dependence, we used Eq. (1) without a thermal- 
reservoir finite-size correction [22, 23]. This correction reduces to replacing T by T — 
Ea/2C in the exponent of Eq. (1), where C is the microcanonical heat capacity of the 
cluster. With C — /cB(3n — 6), where n = 16 is the number of atoms in cubane, the 
closest fit between the modified Arrhenius formula and the numerical-simulation data is 
achieved at Ea — (1.41 ±0.07) eV, which differs substantially from the experimental value 
[8] and, as we will see below, is lower than the height U of the minimum energy barrier 
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to the cubane decomposition. The reason for this difference is unclear, since ear her the 
inclusion of this correction allowed us to describe the experimental data on the fuUerene 
Ceo fragmentation [7]. Unlike fuUerenes C20 and Ceo, cubane consists of unlike atoms 
and the kinetic energy is nonuniformly distributed between the hydrogen and carbon 
subsystems in the thermally insulated cubane [14] . This is hkely to increase the effective 
heat capacity of cubane during its decomposition, and the finite-size correction becomes 
insubstantial. However, this problem requires additional studies. 

Let us find the height U of the minimum energy barrier to cubane decomposition. Figure 
5 presents calculated energies of various isomers of cubane, decomposition products, and 
saddle points on the potential energy hypersurface as functions of the atomic coordinates 
(the details of the calculation procedure can be found in [15, 24, 25]). It is seen that 
the quantity U is determined by the barrier to the transformation of cubane into STCO 
isomer, which completely agrees with the molecular-dynamics simulation data. According 
to our calculations, U = 1.59 eV, which agrees with both the decomposition activation 
energy Ea — (1.9±0.1) eV (which we found from analyzing the numerical simulation 
data) and the experimental value Ea = (1.87 ± 0.04) eV [8]. As we might expect, the 
quantity U is somewhat smaller than Ea, since in experiments (including numerical one) 
cubane can decompose along paths with higher energy barriers. The fact that the barrier 
to the BCT — > COT transition is lower in height than that to decomposition of BCT into 
benzene and acetylene molecules explains why COT is a much more frequent product of 
the cubane decomposition in numerical simulations. 

We also calculated the quantity U by the Hartree-Fock (HF) method without and 
with inclusion of the MoUer-Plesset second-order correction (MP2) and by the density 
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functional method with the B3LYP exchange-correlation functional. All the calculations 
are performed in the 6-31G* basis set. We found that U = 3.95, 3.09, and 3.19 eV, 
respectively. These values are much higher than the experimental value of Ea [8] (despite 
the fact that, as noted above, the inequality U < Ea must take place). Thus, the results 
obtained in the tight-binding model agree much better with the experimental data than 
the first-principle calculations do. This circumstance is due to the fact that we selected the 
model parameters based on the requirement of the best agreement between the theoretical 
and experimental characteristics of the various hydrocarbon molecules [12]. Note that one 
more decisive advantage of the tight-binding model is the fact that the cluster evolution 
over a comparatively long time (on an atomic scale) of i ~ 1 //s can be simulated (whereas 
for the ab initio methods t ^ 1 ps). 

The results obtained in this work have allowed us to find the temperature dependence 
of the lifetime of the metastable cubane CgHg, which may be useful in analyzing possible 
applications of cubane CgRg and sohd cubane s-CsHg as fuel elements. It is also interest- 
ing to develop new methods for synthesizing cubane based on its known decomposition 
products. 
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Fig. 1. Cubane CgHs (schematic). 
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(a) 




Fig. 2. Cubane decomposition products: (a) isomer cyclooctatetraene (COT) and (b) 
CqHq (benzene) and C2H2 (acetylene) molecules. 



12 



(a) 




(b) 




Fig. 3. Isomers of CgHg that form at the cubane evolution stage preceding its decom- 
position: (a) isomer STCO (syn-tricyclooctadiene) and (b) BCT (bicyclooctatriene). 



13 




I IT 

Fig. 4. Logarithm of the hfetime r of cubane CgHg plotted as a function of the reciprocal 
initial temperature T~^: circles are the results of calculations, the solid line is a least- 
squares linear fit, and triangles are the experimental data from [9]. The arrows indicate 
the results obtained in [9] at temperatures T < 474 K and T > 773 K for which, because 
of technical problems, only the lower (r > 40 ms) and upper (r < 0.8 ms) limitations on 
r, respectively, were determined. 
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Fig. 5. Energies calculated by the tight-binding method for various isomers of cubane 
CgHg, its decomposition products, and saddle points Si of the potential energy as a 
function of the atomic coordinates. The cubane energy is taken as a reference point. The 
lines schematically show the paths of possible transitions. The experimental value of the 
cubane decomposition activation energy [8] and the values of the minimum barrier to the 
transformation of cubane into STCO calculated from first principles are also indicated 
(designations of isomers and ab initio methods arc given in text). 
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